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The seniority of debt, which determines the order in which a bankrupt institution repays its debts, 
is an important and sometimes contentious feature of financial crises, yet its impact on system-wide 
stability is not well understood. We capture seniority of debt in a multiplex network, a graph of 
nodes connected by multiple types of edges. Here, an edge between banks denotes a debt contract 
of a certain level of seniority. Next we study cascading default. There exist multiple kinds of 
bankruptcy, indexed by the highest level of seniority at which a bank cannot repay all its debts. 
Self-interested banks would prefer that all their loans be made at the most senior level. However, 
mixing debts of different seniority levels makes the system more stable, in that it shrinks the set of 
network densities for which bankruptcies spread widely. We compute the optimal ratio of senior to 
junior debts, which we call the optimal seniority ratio, for two uncorrelated Erdos-Renyi networks. 
If institutions erode their buffer against insolvency, then this optimal seniority ratio rises; in other 
words, if default thresholds fall, then more loans should be senior. We generalize the analytical 
results to arbitrarily many levels of seniority and to heavy-tailed degree distributions. 

PACS numbers: 89.65.Gh, 89.75.He 


I. INTRODUCTION 

The global financial crisis of 2007-2009 spurred a flurry 
of research on financial contagion using a variety of net¬ 
work models HH23- In most of these models, financial 
institutions interact in just one way, such as making one 
type of loan to one another. In practice, however, banks 
interact in many ways, all of which can contribute to 
a crisis: banks issue loans of different durations and of 
different levels of riskiness; they trade assets with each 
other; and they hold assets in common. Models have only 
recently begun to capture more than one of these kinds of 
interactions [8til0] . Relatively unexplored is the seniority 
of debts, which determines the order in which bankrupt 
institutions repay their debts. Individual creditors pay 
close attention to their assets’ seniority levels, but lit¬ 
tle is known about how the system-wide composition of 
seniority levels affects the risk of large crises. 

A concurrent and mostly independent thread of litera¬ 
ture in physics has studied multiplex networks, or graphs 
with multiple sets of edges (or “network layers”) ; for re¬ 
views, see mm- Variants of percolation on multiplex 
networks have received much attention HM3- Closer to 
financial default contagion are models of threshold cas¬ 
cade on multiplex networks ITMSOj . In these threshold 
models, nodes exist in one of two states (e.g., adopted a 
behavior or not, gone bankrupt or not), and after a node 
adopts the behavior (or goes bankrupt), it influences its 
neighbors in all layers to do so as well. A common 
conclusion from these models of percolation and thresh¬ 
old cascades on multiplex networks is that the system’s 
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vulnerability to large cascades cannot be understood by 
examining the vulnerability of any one layer to large cas¬ 
cades [T5H20] . In particular, in threshold models, the 
interactions among network layers can make large cas¬ 
cades easier (if nodes combine influence from multiple 
types of neighbors in a disjunctive TBJ or weighted [T§] 
way) or more difficult (if nodes combine influence in a 
conjunctive way [20]). 

In this paper, we combine these two strands of research 
to study a financial phenomena of interest: how vulner¬ 
ability to cascading default depends on the seniority of 
debts. Financial institutions (hereafter called “banks” 
for simplicity) hold debt in one another; more generally, 
they have exposure to one another’s default risk through 
interbank lending and/or security holding. Each debt 
contract has a certain level of seniority. When a bank 
goes into default (i.e., when a bank’s liabilities exceed 
its assets), the bank sells its remaining assets to pay off 
whatever liabilities it can in decreasing order of seniority, 
from most senior liability to most junior liability. 

In our model, there exist multiple types of bankruptcy, 
indexed by the highest level of seniority at which a bank 
cannot repay all its debts. Using a technique called the 
cascade condition fTM25] . we study whether a small 
number of initial bankruptcies (of certain types) triggers 
a cascade of many bankruptcies. 

The model has interesting tradeoffs. Self-interested 
banks would prefer to lend at the most senior level in 
order to maximize the assets they would recover in the 
event of its debtors defaulting. However, if all banks 
lend at the most senior level, then the system is espe¬ 
cially vulnerable to large cascades of default, in the sense 
that cascades are likely large for a large interval of net¬ 
work densities. (This extreme case recovers the models 
in HH2].) 
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By contrast, the presence of debts of different seniority 
makes the system less vulnerable to bankruptcies spread¬ 
ing widely. For the case of two seniority levels (called 
“junior” and “senior”), our main result is to derive the 
“optimal” ratio of senior loans to junior loans (for a 
simple class of uncorrelated Erdos-Renyi multiplex net¬ 
works). Here, “optimal” means that it minimizes the size 
of the interval of network densities for which bankrupt¬ 
cies likely spread widely. The optimal multiplex networks 
have on the order of 50% to 100% more senior loans than 
junior ones. As banks’ leverage (i.e., assets divided by net 
worth) increases, this “optimal seniority ratio” increases, 
meaning that more loans should be senior rather than 
junior in order to minimize the size of the set of network 
densities such that cascades are likely large. We also 
show that the optimal seniority ratio is qualitatively sim¬ 
ilar for heavy-tailed degree distributions like those found 
in real interbank lending networks [261 - 128] . These results 
indicate what types of “seniority-multiplex” interbank 
networks are most insusceptible to large cascades of de¬ 
fault. 

The most recent scheme for financial regulation, called 
Basel III, emphasizes the liability structure of individual 
banks, but it does not take into account how seniority 
of debts affect systemic risk. Our results suggest that 
seniority of debt could be an important focus of financial 
regulators and of macro-prudential policies. 

Furthermore, this work contributes to the rapidly 
growing area of research on contagion in multiplex net¬ 
works by introducing a model in which nodes can “ac¬ 
tivate” in more than one way, which we call “multilevel 
contagion”. An interpretation of this model for the 
contagious spread of fads, products, etc. in social net¬ 
works [21] is that the kinds of neighbors a person influ¬ 
ences depends on how excited that person is. Thus, the 
model resembles the multistage model of complex conta¬ 
gion [55] but on multiplex networks. 


A. Related literature 

Theoretical studies of contagion in financial networks 
generally fall into two classes [29]. Some studies con¬ 
sider a particular instance of a financial network, usually 
a complicated and in some cases empirically measured 
object. Other studies consider probability distributions 
over financial networks. 

Our study belongs to the latter class, but two pa¬ 
pers closely related to our analysis belong to the for¬ 
mer m is]- Elsinger ;30] and Gourieroux et al. [.il| 
also assume a senior junior structure of interbank as¬ 
sets. Whereas they obtain a clearing payment vector for 
a given network (as in Eisenberg and Noe [35]), we study 
the cascade condition for ensembles of random networks 
in the same spirit as Watts [21] and followup work in 
physics [TMM [22ti25l 1331135] . The former approach can 
be useful to regulators faced with a real system on the 
verge of a crisis, whereas the latter approach offers ele¬ 
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FIG. 1. (Color online) An example of a bank’s balance sheet. 
Assets (left-hand column) consist of external assets, Is = 4 
many senior interbank loans, and lj = 3 many junior inter¬ 
bank loans. Liabilities (right-hand column) consist of exter¬ 
nal liabilities, bs = 2 many senior interbank liabilities, and 
bj = 3 many junior interbank liabilities. Net worth (a.k.a. 
equity or, for banks, capital) is the difference between assets 
and liabilities. 


gant theory and indicative suggestions |59l. Although fi¬ 
nancial networks change over time, there is evidence that 
their statistical properties (such as degree distributions) 
can be somewhat stable [28] . 


B. Model 

1. Multiplex network of loans of different seniority 

We consider a multiplex network consisting of N nodes 
and M layers. Each node represents a financial institu¬ 
tion (or a “bank” for simplicity). Each layer is a set 
of directed edges (i.e., ordered pairs of nodes) that de¬ 
note loans and, more generally, credit exposures between 
banks. The layer to which an edge belongs is the loan’s 
seniority, and the vector of layers is sorted in increasing 
order of seniority. When a bank goes bankrupt, it re¬ 
pays its liabilities in decreasing order of seniority; that 
is, it repays its most senior liabilities first, its second- 
most senior liabilities second, etc. For the moment, we 
consider a two-layer (“duplex”) network, and the M = 2 
many levels of seniority are called “junior” and “senior” 
(labeled J and S), respectively. 

The number of banks to which a bank lends funds of 
seniority a £ {J, S'} (i.e., the bank’s out-degree in layer 
a) is denoted by l a . Likewise, the number of banks from 
which a bank borrows funds of seniority a (i.e., the bank’s 
in-degree in layer a) is denoted by b a . For simplicity, 
we assume that the sizes of loans (i.e., the notional val¬ 
ues) are the same, so we can consider the networks as 
unweighted. Thus, the out-degree and the in-degree cor¬ 
respond to the bank’s total volume of lending and bor¬ 
rowing, respectively. 

Banks may also hold assets and issue liabilities outside 
the banking system. (For instance, external creditors 
may lend to banks in the form of demand deposits.) The 
balance sheet of a typical bank is illustrated in Fig. [T] 
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2. Multiple types of default 

In a model with two levels of seniority of debt, there 
are two types of defaults (bankruptcies): junior-defaults 
and senior-defaults. A bank is in junior-default if and 
only if its capital (a.k.a. equity or net worth; see Fig. |T]) 
is negative, meaning that it cannot pay all of its junior 
debts. A bank is in senior-default if and only if its capital 
is so far below zero that it cannot pay its all of its senior 
debts nor any of its junior debts. Note that banks in 
senior-default are necessarily in junior-default as well. 

We assume that when a bank goes into junior-default, 
it cannot repay any of its junior-level liabilities, and when 
a bank goes into senior-default, it cannot pay any of its 
junior liabilities nor any of its senior liabilities. Said dif¬ 
ferently, we assume that loss given default of seniority 
a is 100% for interbank liabilities of seniority of a and 
lower. This assumption departs from previous cascade 
models of financial contagion, such as that of Gai and 
Kapadia [T], in which loss given default is 100% for all 
interbank liabilities, no matter their seniority. Table [I] 
makes this assumption concrete: the value of a bank’s 
debt is equal to some positive value (say, 1 unit of cur¬ 
rency) if the bank is solvent (i.e., if its equity is nonneg¬ 
ative); the value of the bank’s junior debt falls to zero 
when the bank becomes insolvent, and the value of its 
senior debt falls to zero when the bank’s equity passes a 
nonpositive threshold value (derived later to be —bj). 


Consider some point in time during a cascade. For a 
certain bank, let mj denote the number of its junior- 
borrowers that have junior-defaulted, and let mg de¬ 
note the number of its senior-borrowers that have senior- 
defaulted. This bank is in junior-default if and only if 
its losses mj + ms exceeds the bank’s equity w, and this 
bank is in senior-default if and only if its losses m j + mg 
exceed the sum of its equity w and its junior liabilities 
bj- 

To exploit the solution technique for threshold cas¬ 
cades in networks P~5H2Ul [22] [23], we express these in¬ 
equalities in terms of losses relative to total out-degree 
h + h'- the bank is in default at layer a £ {J, S} if and 
only if 


mj + m s 

h + Is 


> R a 


Jot if « = J 
ifQ = * 


(i) 


[In Appendix [Aj we generalize the thresholds R a in 
Eq. |l]) to the general case of M > 1 layers.] Figure [2] 
illustrates a bank in junior-default and a bank in senior- 
default in two ways: as losses on a balance sheet (top 
row) and as contagion in a multiplex network (bottom 
row). 

The dynamics of multilevel default contagion can 
be described as a multiplex version of the threshold 
model [21] [2D 33, '36]. (Other multiplex generalizations 
were studied in [T8iH2D] .l The “response functions” for 
each layer are 


TABLE I. Relationship between the equity of a borrower and 
the value of its junior debts and senior debts held by its credi¬ 
tors. Here, we assume that junior and senior debts are equally 
valuable when the borrower is solvent, so values of debts are 
either zero or $1. If the borrower’s equity is negative but suf¬ 
ficiently large (> —bj), then it can pay its senior debts but 
none of its junior debts; the borrower is in “junior-default”. 
For sufficiently negative equity of the borrower, the borrower 
cannot pay neither its junior debts nor its senior debts; the 
borrower is in “senior-default”, and all its debts have value 
zero. 


equity 

status 

junior debt value 

senior debt value 

> o 

solvent 

$1 

$1 

< 0 

junior-default 

0 

$1 

< -bj 

senior-default 

0 

0 


Fj(l,b,m) 


F s (l,b,rh) 


'l if s %ffi L > R J 

0 otherwise 


i if > r j 

0 otherwise 


(2a) 

(2b) 


where l = ( lj 7 lg),b = ( bj,bs ) and to = ( toj,tos )■ If 
the outputs of Fj and Fg are both 0, then the bank is 
solvent; if Fj gives 1, then the bank is in junior-default; 
if Fg gives 1, then the bank is in senior-default (and in 
junior-default). 

We note in passing that a special case of this model 
coincides with a special case of the “multistage complex 
contagion” model of Melnik et al. 25]; for details, see 
Appendix [B] 


3. Thresholds for junior- and senior-default 

Initially, all banks are solvent, and a small fraction 
of banks are chosen to be in junior-default or in senior- 
default. These initial defaults could be caused by the 
loss of external assets. These bank failures cause losses 
for their creditors, who may default as a result, which 
may cause one or more of their creditors to default, and 
so on, resulting in a cascade of default. 


II. MULTILEVEL CASCADES ON 2-LAYER 
MULTIPLEX NETWORKS OF JUNIOR AND 
SENIOR LOANS 

This multiplex model adds a new dimension of com¬ 
plexity to previous threshold models mmmm by 
considering multiple levels of seniority of debt, so to un¬ 
derstand the model we start by simplifying in two other 
dimensions, the distribution of networks and the distri¬ 
bution of thresholds. 
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FIG. 2. (Color online) An example of the balance sheet (and the corresponding multiplex network visualizations) of a bank 
that is (a) in junior-default: mj+ms > Rj(lj + ls) = w and (b) in senior-default: mj + ms > Rs(h+ls) = w + bj. Circular 
and spiky-shaped nodes are solvent and insolvent (a.k.a. bankrupt, in default), respectively, in that layer. Edges point from 
lender to borrower. Senior loans are dashed green arrows; junior loans are solid red arrows. In panel (a), the bank labeled A 
suffers losses of mj + ms = 2 on its assets side (assets are listed on the left-hand side of the bank’s balance sheet) because 
mj = 1 of its (lj = 3 many) junior debtors are in junior-default and because ms = 1 of its (Is = 4 many) senior debtors are in 
senior-default. These losses exceed bank A’s equity ( w ), so bank A cannot pay off all its (bj many) junior liabilities. Because of 
the assumption that loss given default at layer a is 100% for all creditors in that layer, all of A’s junior creditors do not receive 
repayments on their loan to A, and we say that A goes into junior-default. However, bank A does not go into senior-default 
because it can still pay its senior liabilities (i.e., because mj + ms < w + bj). In panel (b), bank A suffers so many losses 
(mj + ms = 5) that it cannot repay its junior creditors nor its senior creditors, so bank A goes into senior-default. 


First, each layer is modeled as an independent, lo¬ 
cally treelike, random directed network generated by the 
configuration model. The sparsity and independence 
of the layers implies that the number of “overlapping 
edges” (i.e., junior and senior loans between the same 
pair of banks) is negligibly small for sufficiently large net¬ 
works [57] , 

Second, following previous studies m , we assume 
that every bank has the same junior-default threshold 
Rj = w/(lj + Is)- This assumption implies that all 
banks have the same ratio of capital to total interbank 
assets. In the real world, banks’ capital-to-asset ra¬ 
tios are not identical, of course, but they are somewhat 
similar [38]. It is straightforward to consider hetero¬ 
geneous capital-to-asset ratios by considering heteroge¬ 
neous thresholds [3] [22] ■ 


A. Analytical approximation 

Let 4>q and (f >g denote the probabilities that a node is 
initially in junior-default or in senior-default, where 0 < 
cj) j[<Cl and 0 [0, <^g]. Next we adapt to this model a 
technique for calculating whether cascades likely spread 
widely, the cascade condition , introduced in and 

generalized to other multiplex network models in T5U2U] . 


1. Recursion equations 

Consider a bank chosen uniformly at random, and let 
the network hang down like a tree from this root node 
by doing a breadth-first search along the directions of the 
edges. We define as the probability that a node t +1 
hops from the bottom of the tree is in a-default because 
sufficiently many of its debtors (i.e., its out-neighbors) 
are in default. The nodes at the bottom of the tree are 
initially in a-default with probability </>q. To compute 
from </>“, we condition on whether a node located t+ 
1 hops above the leaves of the tree is initially in a-default, 
which occurs independently with probability <f> q. If this 
node is not initially in a-default, then we condition on 
the number of debtors of this node that have defaulted. 

For a sufficiently large graph, the fractions of nodes 
in junior- and senior-default at the end of the cascade 
are well approximated by a fixed point </>^) of the 
recursion equations 

ff+i =9 (a \ri,tf) 

= ft+ (1-4%) £ E^pf/o^o^ow 

1 bj 

h Is 

X E E (3) 

mj—0 ms—0 
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r _ r t o'! i ck, borrow j a,loan j , ,1 

for a G {J, o}, where and p L denote the 

in- and out-degree distributions, respectively, on layer 
a £ {J, S}. In Eq. we approximated the numbers 
of neighbors in junior and senior default as independent 
binomial random variables, the probability mass func¬ 
tions of which are written using the shorthand notation 
B l m (() f>) = 1 — <j>) l ~ m \- This approximation works 

well when the network layers are independent and hence 
the number of overlapping edges is negligibly small E2- 
(If overlap is significant, then one could adapt the tech¬ 
niques introduced in a study of mutual percolation [TB] 
to this threshold model.) Variants of Eq. [3] appear in 
several studies of this threshold model; see, for exam¬ 
ple, [221 Eqs. (1)—(3)], [231 Eqs. (1) and (4)], [El Eqs. 
(1)—(3)], [El Eqs. (11)—(13)], |251 Eqs. (A9)-(A11)], [El 
Eqs. (1)—(3)] and, for a rigorous treatment, [7] Theorem 
1]; the closest variant is the study of directed (but not 
multiplex) networks [231 Eqs. (14) and (15)]. 

In short, we approximate the multiplex graph as a tree 
according to outgoing edges, and is the probability 
that a bank (located t+1 hops from the leaves of the tree) 
goes into a-default because it was initially in a-default or 
because of the defaults of its debtors. Note from Eq. © 
that going into senior-default due to defaulted debtors 
implies having gone into junior-default due to defaulted 
debtors [i.e., Fs(l,b,fh) = 1 implies Fj(l,b,m ) = 1], so 
the fraction in senior-default never exceeds the fraction 
in junior-default [i.e., </>f < <j>{ for all t > 0]. 

For convenience, we rewrite the recursions ([3]) as 





(4) 


Iterating Eq. Q to a fixed point (</>£,, ^) gives the ex¬ 
pected fractions of banks in junior-default and senior- 
default, respectively, at the end of the cascade. 


2. Multiplex cascade condition 


Now we obtain a simple expression that approximately 
captures, in the limit of network size approaching infinity, 
for which parameters a vanishingly small seed triggers a 
cascade that results in a finite fraction of banks in junior- 
or senior-default. The first-order cascade condition is the 
linear instability of Eq. Q at the origin (0,0), which 
provides a sufficient condition for the numbers of junior- 
and senior-defaulted banks to grow. The Jacobian matrix 
of the right-hand side of Eq. 0 is denoted by 

J (5) 


The first-order cascade condition is that the largest eigen¬ 
value of the Jacobian matrix, A max (j7’), exceeds 1: 


A 


max 


(J) 


Jjj + Jss + \J ( Jjj - Jss) 2 + 4 JjsJsj 


> 1 . 


( 6 ) 


In Appendix 0 we simplify the Jacobian matrix and 
show that, for this model, its largest eigenvalue equals 
its trace. 


B. Cascade regions 

If many banks are in junior-default, the fraction of 
banks in senior-default could be large or small. Thus, 
there are two other cascade conditions: one for junior- 
defaults and another for senior-defaults. Next we derive 
these conditions in some generality (for ensembles of mul¬ 
tiplex networks with independent layers that are speci¬ 
fied by their degree distributions). We also illustrate all 
three cascade conditions (junior, senior, and multiplex) 
for a particular case of interest: a multiplex network of 
two uncorrelated Erdos-Renyi networks of mean degree 
(lj) and (Is), respectively. 


1. Cascades of junior-default (ignoring senior-default) 

First consider contagion only on junior loans. The 
“junior-only” cascade condition is that a junior-default 
causes on average at least one more junior-default, or 
equivalently that the top-left entry of the Jacobian ma¬ 
trix, Jjj , exceeds 1. In Appendix [C] we show that this 
junior-only cascade condition is the inequality 

Jjj =E[Ij(1 1>Rj ( 1j+ I s ))] > 1. (7) 

The set of parameters satisfying this junior-only cas¬ 
cade condition 0 is shown in orange in Fig. [3](a) for 
uncorrelated Erdos-Renyi networks of mean degree (lj) 
and (Is)- If the senior-layer is the empty graph (i.e., if 
(Is) = 0), then we recover standard cascade models on a 
single-layer network HH31 HQ 123 \M- Creating more se¬ 
nior loans (i.e., increasing (Is)) shrinks the “junior-only” 
cascade region. This shrinking occurs because a bank’s 
equity w = Rj(lj+ls) increases with its number of senior 
loans Is [recall Eq. Q]. Thus, as the expected number 
(Is) of senior loans per bank increases, a typical bank’s 
buffer against junior-default (w) increases, which makes 
a junior-default less likely to cause other junior-defaults. 


2. Cascades of senior-default (ignoring junior-default) 

In a similar way, consider contagion along senior loans 
only. The “senior-only” cascade condition is that a 
senior-default causes on average at least one more senior- 
default, or equivalently that the bottom-right entry of 
the Jacobian matrix, Jss, exceeds 1. In Appendix[Cj we 
show that this senior-only cascade condition is 

JSS = E[ls(^l-bj>Rj(lj+ls))] 

= Po’ borrow E[Zsl (8) 


2 
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The first-order cascade condition considers defaults trig¬ 
gered by the default of just one out-neighbor (i.e., de¬ 
faults due to mj + m-s = 1); thus, it captures senior- 
defaults of banks that lack junior debts (bj = 0), which 
explains the factor pA borrow j n Eq. (jsj) . We found that 
capturing defaults caused by two bankrupt creditors (us¬ 
ing the second-order cascade condition) was not neces¬ 
sary to approximate simulations well, as shown in Ap¬ 
pendix [D] and in Fig. [bj 

The parameters satisfying this senior-only cascade con¬ 
dition ^ are shown in green in Fig.[3])a) for uncorrelated 
Erdos-Renyi networks of mean degree (lj) and (Is)- This 
region shrinks with the density of junior loans (lj) for the 
same reason that the junior-only cascade region shrinks 
with (Is)- 

3. Multiplex cascades 

The blue region in Fig.[3](a) is the multiplex cascade re¬ 
gion, the set of average out-degrees ((lj), (Is)) such that 
the first-order cascade condition [inequality |g])] is satis¬ 
fied. Note that, as pointed out in IH page 3], we have 
A ma x(i7) > max{ Jjj,Jss}- Thus, a sufficient condition 
for the multiplex cascade condition to hold is that at least 
one kind of default a € {J, S'} is supercritical (meaning 
that J aOL > 1) even in the absence of the other kind of 
default. Thus, the blue region in Fig.^a) contains both 
the orange and green regions. 

Importantly, however, the blue region contains pairs of 
intermediate edge densities ((lj), (Is)) outside the junior- 
only and senior-only cascade regions. In this part of pa¬ 
rameter space, neither junior-defaults nor senior-defaults 
alone suffice to cause global cascades (i.e., Jjj < 1 and 
Jss < 1), but a large cascade of defaults can neverthe¬ 
less occur. These cascades, akin to the “multiplexity- 
facilitated cascades” in [18] [20], result from the interac¬ 
tion of junior- and senior-defaults. 


C. Optimal seniority ratio 

Note that the multiplex cascade region, show in blue 
in Fig.^a), is not symmetric with respect to the average 
out-degrees in the two layers, (lj) and (Is)- This asym¬ 
metry of the cascade region results from the asymmetric 
role that junior and senior loans play in the response 
function [Eq. §]■ Given this asymmetry, what ratio of 
junior and senior loans is optimal in some sense? 

Suppose that the average total degree ( lj) + (Is) is 
kept constant and that we tune the fraction of loans 
that are junior. For some loan densities (lj) + (Is)-, 
the multiplex cascade condition © holds no matter the 
fraction of junior loans. For example, the line segment 
{((lj),(ls)) ■ 4 = (lj) + (Is)} is contained within the 
multiplex cascade region in Fig. [3](a), so defaults spread 
widely no matter the composition of seniority levels. In 
other cases, such as the line segment of mean degrees 



expected number of junior loans per bank = (lj) 



seniority ratio ( l s ) / (lj) 

FIG. 3. (Color online) (a) The multiplex cascade region 
[the edge-density parameters ((h), (Is)) satisfying Eq. jfjj, 
drawn as a blue region with a thick, solid boundary] con¬ 
tains the junior-only and senior-only cascade conditions given 
by Eqs. 0 and |8| and drawn as an orange region (with a 
dashed boundary) and green region (with a dotted boundary), 
respectively. Moreover, the multiplex cascade region is asym¬ 
metric in the densities of junior and senior loans, ((h), (Is))- 
The optimal seniority ratio, defined in Eq. |9|, is drawn as a 
black line. This ratio (Is)/(h) ^ 1-79 minimizes the cascade 
window, the set of mean degrees satisfying the first-order cas¬ 
cade condition, Eq. ©■ The case of equally many junior and 
senior loans (on average) is shown with a dashed line, which 
has a slightly larger cascade window, (b) The cascade win¬ 
dow as a function of the ratio of senior loans to junior loans, 
(Is)/(h), achieves a minimum size at « 1.79. 


satisfying 7 = (lj) + (Is) in Fig.[3](a), the cascade condi¬ 
tion © holds if and only if the fraction of junior loans 
is very small or somewhat large. In cases like this one, 
having an intermediate fraction of junior debts prevents 
defaults from spreading widely. 

Now suppose that we can set the ratio (Is)/(h) of 
the average numbers of senior and junior loans per bank, 
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but the total density of loans (Is) + (Is) can vary. For 
instance, the number of loans may increase during an 
economic boom and decrease during a recession, yet the 
fraction of loans that are senior does not change. Thus, 
we are considering rays through the origin in Fig. [3](a). 

For “seniority ratios” a > 0, define the cascade win¬ 
dow to be the set of average degrees ((Is), (h)) such that 
(Is)/(h) = o and such that the first-order cascade condi¬ 
tion 0 holds. Figure^b) shows this cascade window as 
a function of the seniority ratio (Is)/(h)- Note that the 
size of the cascade window is minimized at a unique value 
of (Is)/(h), which is approximately 1.79 in Fig. [3j We 
define the optimal seniority ratio to be this minimizer, 
or more precisely, 


M 

(h) 


optimal := argmin |{r > 0 : A max 

cr>0 


(J)> i}| 


(9a) 


where r is the distance from the origin, | • | denotes 
Lebesgue measure on R, and A max (i7) is evaluated at 


(h) = 



(h) 


ra 

\/l + a 2 ’ 


(9b) 


where we have changed from polar to Cartesian coor¬ 
dinates via (r, arctan ct) i —> ((h), (h))- If this optimal 
seniority ratio ((Is)/(h)) optimal is greater than one (as 
in Fig. [3]) , then senior debts should be issued more than 
junior debts (on average) for the sake of minimizing the 
size of the set of network densities such that bankrupt¬ 
cies spread widely. That is, if a financial network keeps 
approximately the same ratio of senior to junior loans 
((Is)/(h)) hut changes the total number of loans (say, 
over the course of a business cycle), then there should be 
more senior loans than junior loans in order to minimize 
the risk of a large cascade of defaults. 

Intuitively, senior loans are stabilizing because they are 
less significant early in a cascade: a node in junior-default 
causes losses (of assets) for its lenders (in-neighbors) in 
the junior layer only. Thus, changing some loans from ju¬ 
nior to senior makes large cascades less likely, and hence 
it shrinks the cascade window. However, not all loans 
should be senior, because that case is identical to all 
loans being junior, which has a dangerously large cas¬ 
cade window [look at the horizontal and vertical axes of 
Fig. l: a)]. In between those extremes lies an optimal ra¬ 
tio of senior to junior loans. That optimal ratio exceeds 
one [the dashed black 45-degree line in Fig. |3](a)] because 
senior loans are less contagious, in the sense that nodes 
in junior-default can still repay their senior debts. 

In Fig. [7] of Appendix [DJ we show that the three 
types of cascade regions shown in Fig. [3ja) (junior-only, 
senior-only, and multiplex) closely agree with the fixed 
point of the recursion equations @ and with numerical 
simulations of the model for uncorrelated Erdos-Renyi 
multiplex networks. These comparisons also show that 
within the multiplex cascade region many nodes end up 
in junior-default, but many nodes are in senior-default 
only when senior loans greatly outnumber junior loans. 


In practice, the number of defaults of any type is of pri¬ 
mary concern, so our definition of the optimal seniority 
ratio in Eq. @ ignores types of defaults. 


D. Effect of varying the thresholds (ratio of equity 
to interbank assets) 

How does the cascade region (and its optimal senior¬ 
ity ratio) change as the banks’ threshold Rj changes? 
Here we keep the assumption that banks all have the 
same junior-default threshold Rj , but we vary Rj. The 
multiplex network is still a pair of uncorrelated Erdos- 
Renyi random graphs with mean degrees (lj) and (Is)- 
If banks have a smaller capital ratio Rj , then default be¬ 
comes more likely, so the multiplex cascade region grows, 
as shown in Fig. |4|a). Conversely, as Rj increases, the 
cascade region shrinks, and it can even splinter into two 
regions [see the purple regions for Rj = 0.25 in Fig.[4|a)]. 
In cases like this one with two disjoint cascade regions, 
the optimal ratio defined in Eq. § does not exist because 
the set {r > 0 : A max (J') > 1} is empty for some a. [A 
sensible definition for the optimal seniority ratio could 
be the slope that maximizes the distance from the ray 
through the origin (with that slope) to the two cascade 
regions.] 

The asymmetry between junior and senior loans also 
changes with the junior-default threshold Rj. In particu¬ 
lar, the optimal seniority ratio increases as Rj decreases, 
as shown in Fig. |4][b). 

In some financial crises, many banks find themselves 
in financial distress, because of losses in external assets, 
for instance. The model studied here approximates this 
scenario as a small junior-default threshold Rj for every 
bank (and a small fraction of banks initially in default). 
Decreasing Rj enlarges the cascade region, as illustrated 
in Fig. |4][a). However, the banks can reduce systemic 
risk if they increase, on average, the fraction of debts 
that are senior. Figure |4|b) shows the optimal amount 
of that increase. 


III. MORE THAN TWO SENIORITY LEVELS 


Now we show that mixing among several seniority 
levels can further reduce systemic risk (in the sense of 
shrinking the cascade region). 

We generalize the model to M > 1 seniority levels. 
Let (1,2,..., M ) denote the indices of the seniority lev¬ 
els in increasing order of seniority, and let I?i denote 
the threshold for default at the most junior layer. The 
response function for default at level i, derived in Ap¬ 
pendix [A] is 


Fi(l, b, m) 


1 if Y,?=irn,-TirJibk>Rii:? = 1 l, ' 

0 otherwise 

( 10 ) 







(a) 14 


first-order cascade region 
■ */ = 0.1 

■ Rj — 0.15 

■ R, = 0.2 

■ Rj = 0.25 


2 4 6 8 10 12 14 

expected number of junior loans per bank = (lj) 



junior-default threshold Rj 

FIG. 4. (Color online) If banks hold less capital (equity) w 
relative to their interbank assets lj + Is (he., if the threshold 
for junior-default, Rj, decreases), then (a) the cascade re¬ 
gion (the set of parameters satisfying the first-order cascade 
condition Q] enlarges, and (b) the optimal ratio of senior to 
junior loans increases. The first-order cascade condition |6| 
considers defaults that occur due to one defaulted neighbor 
(for details, see Appendix [Cf. Thus, when Rj falls below the 
inverse of an integer, 1 /k, nodes with total degree k and with 
one defaulted neighbor will junior-default, so the cascade re¬ 
gion changes at Rj = 1/fc, and hence the optimal ratio jumps 
up at Rj = 1/fc. 


By convention, sums of the form Y^t =i are zer0 - 


A. First-order cascade condition 


ceeding one. The (*, j) entry of is 


-E m h 1 i-j: i k - 1 1 b k >R 1 (j:™ =1 l a ) 


^k ,borrow 


M il a ) 


i -1 

=n^o 

k =1 

i -1 M 

n k, borrow \ 1 T s, loan 7 

Po 11 Pis ?■ 

fc=1 i : E"=ih<i/iii s=1 


(11a) 

(lib) 

(11c) 


Here, products of the form J|° =1 are defined to be 1, 

• ttp l.loan M, loan i 

and the expectations Em are over p (| ,... and 

1,borrow M— 1,borrow mi • • fTn" \ 

over p bi ,... ,Pb M _ 1 • ine expression m lib) 


equals that in (11a) because the threshold R\ > 0, so 


the indicator in (11a) is 1 only if the bank has no liabil¬ 
ities of seniority k (i.e., bk = 0) for all k = 1,2,...,* — 
1; furthermore, because a bank’s in- an d out -degrees 
are independent, the indicator in Eq. (11a) factors 
to lh )b =ovfce{i, 2 ,...,i-i}li>_R 1 (^y_ l i a )- Equation fllTc] ) 


writes out the expectation in Eq. (lib) by summing over 


Z e {0,1,2,... } M such that J2a=i < V-Ri- 

Observe from Eq. (lib) that J M consists of rows that 
are linearly dependent. Thus, J M has rank one and only 
one non-zero eigenvalue. Because the trace of a matrix 
equals the sum of its eigenvalues, t he lar gest eigenvalue 
of J M is its trace, which, from Eq. (11c), is given by 


M i— 1 


M 


tr J M = Y, n %’ borrOW 


e n 


Pi 


s,loan 


k. (12) 


i—1 k—1 


l:E"ih<l/iil ' 


Therefore, we can analytically express A max (J M ) by us¬ 
ing the degree distributions and the threshold value Ri 
on the most junior layer. This expression immediately 
gives the first-order cascade condition, trJ M > 1. 


B. Distinguishing more seniority levels reduces risk 
of large cascades of default 


If we realize that not all loans have the same se¬ 
niority level, how does vulnerability to large cascades 
change? For simplicity, we consider “splitting” an Erdos- 
Renyi random graph into M independent, identically dis¬ 
tributed layers (as in DSD- That is, we compare a single¬ 
layer Erdos-Renyi random graph with mean out-degree 
z to a multiplex Erdos-Renyi random graph with mean 
out-degree z/M in each of the M layers. 

The Jacobian matrix J’ AI ~ EK for an M -layer directed 
Erdos-Renyi random graph with edge density (li) in each 
layer i £ {1, 2,..., M} is, from Eq. (11), 


The cascade condition, a straightforward generaliza¬ 
tion of the one for M = 2 layers in Sec. |II A 2[ is that the 
Jacobian matrix J M £ I MxM has largest eigenvalue ex- 


stM-ER 


r(rWI -i) 


exp 


' 3-1 

-E>> 

. k =1 


(13) 
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where the incomplete gamma function T (x, y ) = 

u x ~ 1 e~ u du and the gamma function r(a:) = F(x, 0 ). 

Figure [5] shows the resulting cascade region (i.e., the 
parameters satisfying trj7 M ' ER > 1) for M e {1,2,3,4} 
many Erdos-Renyi layers that each have mean out-degree 
(li) = z/M. Here, the parameters are the threshold R\ 
for default at the most junior level (horizontal axis) and 
the mean total out-degree 0 (vertical axis). 


number of seniority levels (M) grows. It is possible that if 
there were more seniority levels, then banks would erode 
their buffers against insolvency, R\ , thereby pushing the 
system back toward the cascade region, where bankrupt¬ 
cies likely spread widely. 


IV. HEAVY-TAILED DEGREE 
DISTRIBUTIONS 


first-order cascade region 
M = 1 seniority level 
M = 2 seniority levels 
|M = 3 seniority levels 
■ M = 4 seniority levels 



0.10 0.15 0.20 0.25 0.30 

threshold R t in the most junior layer 


0.35 


FIG. 5. Distinguishing more types of seniority reduces the 
size of the cascade region. In other words, with more types of 
seniority, the expected size of a default cascade is negligibly 
small for more pairs of thresholds Ri and expected number 
of loans per bank z = Y/iLi^i)- Here we consider splitting a 
directed Erdos-Renyi network into M £ {1,2, 3, 4} many in¬ 
dependent, identically distributed layers, each with the same 
mean degree (k) = z/M. 


Notice in Fig. [5] that the cascade region shrinks with 
the number M of seniority levels, meaning that the sys¬ 
tem is less vulnerable to large cascades of default. (By 
contrast, in the model in [18j . the cascade region grows 
with M.) For some capital-to-interbank-asset ratios R\, 
the cascade region is even eliminated as M grows (e.g., 
Rx = 0.3 for M = 1). 

The reason the cascade region shrinks with M is that 
making a loan more senior makes it less significant early 
in a cascade. For instance, if M = 1, then a defaulted 
node reduces the assets of all of its lenders (in-neighbors), 
whereas if M >1, then only the in-neighbors in layer 1 
lose an asset, so less contagion occurs. The decline in 
importance of more senior layers is captured by the mul¬ 
tiplicative factor nr=U’ b ° rrOW i n Eq. ( 11 c) and by the 
exponential factor in Eq. ©■ In short, distinguishing 
more types of seniority makes more interbank debts in¬ 
significant (early in a crisis); hence a systemic crisis is 
less likely. 

This result hinges on the assumption that banks’ 
buffers against insolvency R± remain constant as the 


Real networks of loans among banks show evidence 
of somewhat heavy-tailed degree distributions [ 261128 ] . 
For instance, the network of Brazilian banks in 2007 
2008 [28] has in- and out-degree distributions that have 
tails approximated by a power law oc fc -7 with exponent 
7 ss 2 . 8 , and the mean in- and out-degree are both ap¬ 
proximately 8.5. 

To capture such heavy-tailed degree distributions, we 
use the static model of scale-free networks GUIS)], which 
allows tuning the mean degree and the exponent 7 of the 
degree distribution [44]. Because we lack access to data 
on seniority levels of loans, we first create a random graph 
with a certain mean out-degree, and then we assign a 
fraction of directed edges (chosen uniformly at random) 
to be junior, and the remaining are senior. 

Figure [ 6 ] shows that, in numerical simulations and in 
the analytical approximation ([ 6 ]), the fraction of loans 
that are junior can determine whether large cascades of 
defaults tend to occur. Notice the “valley” in Fig. [ 6 ] 
there is a narrow range of “optimal” values of the frac¬ 
tion of junior loans (horizontal axis of Fig. 6 | that shrinks 
the set of mean-degrees such that cascade size is large. 
This optimal fraction of junior loans is approximately 
0.30, which translates to an optimal seniority ratio of 
(Is)/(h) ~ (1 — 0.30)/0.30 « 2.3, which is qualitatively 
similar to the optimal seniority ratio 1.79 found for mul¬ 
tiplex Erdos-Renyi networks in Fig. [3] 

Notice in Fig. [ 6 ] that the empirical mean out-degree in 
the Brazilian network (approximately 8.5) [25] intersects 
the black regions (where the mean cascade size is large) 
whenever the fraction of junior loans is large or small. We 
emphasize that we do not know the capital levels of these 
Brazilian banks, which affects the threshold for junior- 
default, Rj. However, Fig.[ 6 ]suggests that seniority levels 
of loans in real-world networks could potentially affect 
whether default might spread widely or not. 


V. CONCLUSION AND DISCUSSION 

We model seniority of debt using a multiplex network, 
and we study the likelihood of bankruptcy spreading 
widely using a generalized cascade condition. The model 
has interesting tradeoffs. Whereas self-interested banks 
would prefer to make their loans at the most senior level, 
a mixture of junior and senior debts optimally reduces 
systemic risk, in that it minimizes the size of the set 
of network densities such that cascades are likely large. 
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on o? n.4 or os in 



fraction of loans that are junior 


FIG. 6. (Color online) For heavy-tailed degree dis¬ 

tributions like those in the Brazilian interbank lending 
network [28], having « 30% of loans be junior (and the rest 
senior) reduces systemic risk, in that it minimizes the size 
of the set of mean total out-degrees such that bankruptcy 
spreads widely. The gray-scale background shows the 
average fraction of banks in junior-default (top panel) and in 
senior-default (bottom panel) in simulations on 2400 banks, 
begun from 10 banks initially in senior-default, averaged 
over 30 simulations. The green lines mark the boundary of 
the casca de region [A max (J) = 1 with J computed from 
Eq. (C41]. The junior-default threshold of every bank is 


Rj = 0.18, as in Fig. [3 
the static model no 


Directed graphs are generated with 
with out-degree distributions having 
tail behavior oc k 1 with 7 = 2.83, which approximates the 
value found in the Brazilian interbank lending network [28] . 
A fraction of loans (horizontal axis) chosen uniformly at 
random are junior, and the rest are senior. The “optimal” 
fraction of junior loans (that minimizes the height of the 
green line) is ss 0.30. For the mean out-degree (vertical 
axis) of £3 8.5 observed for Brazilian banks [23], the expected 
cascade size is small if and only if an intermediate fraction of 
loans are junior. 


away, such as interest rates for different seniority levels, 
which affect the incentives of banks. 

Traditionally, cascades of losses in financial networks 
have been studied using the Eisenberg-Noe algorithm [32] 
for computing a market-clearing payment vector (with 
losses shared proportionally), which was recently ex¬ 
tended to multiple seniority levels [51] . The Eisenberg- 
Noe algorithm is economically meaningful in normal 
times, but as Cont et al. |2B] point out, the algorithm is 
not so reasonable for modeling an extreme event such as 
the default of a large bank. We instead take an approach 
typical in physics of studying ensembles of financial net¬ 
works mm- We show that cascades on networks with 
different seniority may be analyzed with a variant of mul¬ 
tiplex threshold cascade models studied in physics pa- 
120 ]. 

The seniority of debt is just one dimension of multi- 
plexity in financial networks. Other dimensions of debt 
include the maturity of loans and whether loans are col¬ 
lateralized M- Furthermore, other important interac¬ 
tions include the trade of financial assets (43j and com¬ 
mon exposures (i.e., multiple institutions holding the 
same financial assets) [HI SI] ■ 

We call our model “multilevel contagion” because a 
node can become “contagious” in one layer without nec¬ 
essarily becoming contagious on another layer (e.g., a 
bank can be in junior-default but not in senior-default). 
By contrast, in other models of threshold contagion in 
multiplex networks [18M20] and in variants of percola¬ 
tion on multiplex networks HMZ], when a node changes 
state in one layer, it necessarily changes state in other 
layers. Multilevel contagion might give insight into other 
kinds of social contagion. For example, one might use it 
to model the adoption of a new product if “fans” of the 
product spread influence to neighbors of one type but 
“super-fans” spread influence to neighbors of two types 
(a multiplex version of the multistage complex contagion 
model [25]). We hope our model will stimulate further 
research on these many types of cascades in social, finan¬ 
cial, and other multiplex networks. 


Such disparity between individual interests and system- 
wide stability is the “regulator’s dilemma” [42]. For two 
types of seniority, and for edges placed uniformly at ran¬ 
dom (i.e., uncorrelated Erdos-Renyi networks) or accord¬ 
ing to a heavy-tailed degree distribution, we show that 
the optimal multiplex networks have on the order of 50% 
to 100% more senior debts than junior debts. To give 
theoretical support for policy suggestions, however, will 
require considering other factors that we have abstracted 
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Appendix A: Thresholds of default for M levels of seniority 

Here we consider the general case of M > 1 levels of seniority of interbank debt (rather than just two, junior and 
senior), and we label the levels of seniority by a € {1,2 ,,M}. We choose the convention that debts of seniority a 
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are more senior than debts of seniority a' whenever a > a'. 

Recall from Fig. [T] the stylized balance sheet of a bank: a bank’s assets consist of external assets e and l a many 
interbank loans of seniority a £ {1,2,..., M}, while its liabilities consist of external liabilities d and b a many interbank 
liabilities of seniority a £ {1, 2,, M}. We assume that the external liabilities d (which could include, for example, 
demand deposits) are more senior than debts to other financial institutions. The equity w is the difference between 
assets and liabilities, 


M 

w = e~d+ "^2(l a - K)- 

a=l 

Now suppose that this bank has some interbank assets (i.e., la > 0) and that it suffers losses on m a of its 

l a many interbank assets (for a £ {1, 2,..., M}), so that its interbank assets now have total value J2a= i(fc — m a ). 
Then exactly one of the following events occurs: 

• The bank is solvent if and only if the capital buffer w can absorb the losses, i.e., w > Y2!= i m «- 

• The bank is in default at layer a! £ {1,2 ,,M} but not in any layer > a' + 1 if and only if 

a.'— 1 M ol 

w + E b a < E m a < w + E b a . (Al) 

a—1 a—1 a=1 

• The bank is in complete default , meaning that it cannot pay off any of its interbank debts nor any of its liabilities 
to external creditors, if and only if 


M 


M 


^ m a > w + ^ b a . 


Ot — 1 


a—1 


Here, we focus on cascades in the interbank network; for this purpose, we do not need to distinguish between default 
at the most senior level (M) from complete default because both have the same effect on the banks that had lent to 
the defaulted bank. 


Rearranging the strict inequality in Eq. (Al) and dividing by the total lending XEi la > 0 gives the response 
function in Eq. (10). 


Appendix B: Relationship with “multistage complex contagion” 

Our model is related to the “multistage complex contagion” model of Melnik et al. [25] . In their model, nodes 
exist in one of three states: inactive, active, and hyper-active. Their hyper-active nodes are active nodes with greater 
influence on their neighbors. To see the commonality between the models, consider our model with every node 
having the same junior threshold Rj and with the multiplex network consisting of two identical layers (i.e., all the 
multiplex edges are “fully overlapped” [1611371 ). Then every neighbor in junior-default (respectively, in senior-default) 
contributes mj + ms = 1 (respectively, mj + ms = 2) in the numerator of the response function. Thus, this model 
is equivalent to a special case of the model in [25| in which i) the “bonus influence parameter” /3 (which captures the 
extra influence of hyper-active nodes) takes the value 1; ii) the network is a directed single-layer graph; and iii) the 
thresholds for activation and hyper-activation are Rj and Rj + bj/[lj + Is) = Rj + bj/(2lj), respectively. To model 
real financial systems, however, edges should be only partially overlapped, a phenomenon studied in the contexts of 
statistical mechanics in m and of percolation in mi- Our multilevel contagion model could be generalized to such 
a partially overlapped case by adapting techniques from the study of mutual percolation on multiplex networks with 
arbitrary amounts of overlap m, a nontrivial task left for future work. 


Appendix C: Derivation of the first-order cascade condition 

Now we obtain a simple expression that approximately captures for which parameters a vanishingly small seed leads 
to a cascade that results in a finite fraction of banks in junior- or senior-default. Note that if the threshold Rj > 0, 
then (0, 0) is a fixed point of the recursion Q [because B^( 0) = l m= o and F a (l, b, 0) = 0 provided that Rj >0]. The 
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first-order cascade condition is the linear instability of Eq. Q at this fixed point at the origin ((/>/, (ft) = (0, 0) in the 
double limit (f$ —> 0 and </g —» 0. This linear instability provides a sufficient condition for the numbers of junior- and 
senior-defaulted banks to grow (at least initially). The Jacobian matrix of the right-hand side of Eq. 0 is 


f Jjj Jjs\ 

\Jsj Jss) 

\ E [^j(ll-bj>,Rj(Zj+Z s ) &j>i?j(Zj+Zs))] E /s(-®-l — bj>Rj(l j+ls) ^ — b j>Rj(l J+ls) )] 


(Cl) 


(C2) 


(The expectations E are over the degree distributions ° an ,p /g ° al \Pbj ° 1Iow .) Assume that Rj > 0, so that the 
factors lojjj in the top row of (C2) vanish. Also, note that bj > 0 and that Rj > 0 and lj + Is > 0, so the 
factors 1 -bj>R,,(ij+i s ) in the bottom row of Eq. (C2) vanish. Finally, the other indicator in the bottom row of 
Eq. (C2), namely simplifies to lbj=oli>fij(jj+z s ) because bj £ {0,1,2,...}. Because each node’s 

in- and out-degrees ( lj,bj,ls,bs ) are all independent, we can factor the expectation E [^alb J =gl 1> _R J (; J +i s )] = 
Pg’ b °rio w E for both a £ {J,S} (he., for the bottom-left and bottom-right entries of Eq. (C2|). It 

follows that 


rr _ ( E [^!l>Jij(Jj + is)] E ^s!l>.Rj(Zj + Zs)] \ 

J -{p J 0 ’ h —E[ljt 1>RAlj+ls) ) p J 0 ^E[lsl 1> R Alj+ls) }J 

= E [^jll>iJj(Zj+Z s )] borrow + E [^sll>i?j(Zj+Z s )] ^J ,borrow 

Note that det J = 0. Thus, from Eq. ([o]), the largest eigenvalue of the Jacobian is A max (i7) = tr J = Jjj + Jss- 
It should be pointed out that this equality, A max ( l 7) = trJT, does not necessarily hold true in other variants of 
multiplex contagion models. For example, in multiplex models with (independent) undirected networks and/or an 
exogenous senior-default threshold Rs, one would generally have A max (J) ^ tr J. A more detailed explanation is 
given in Appendix |E| 


(C3) 

(C4) 


Appendix D: Comparison between the first-order cascade condition with the fixed point of the recursion 

equations and with numerical simulations 

Figure [7] shows that the first-order cascade condition [the parameters satisfying Eq. ([6]), which are enclosed by the 
solid blue line in Fig. [ 7 ] agrees closely with the fixed point (</>/[,, ) of the recursion equations @ [Fig- [0(a) and^b)] 

and agrees reasonably closely with the results of numerical simulations [Fig. 0c) and [70)]. In the recursion equations 
and simulations, we use the same initial condition, with </>g = 4>q = 5 x 10 -4 . The threshold for junior-default is 
Rj = 0.18, as in Fig. [3j 

Notice in Fig. [T] that when a large cascade occurs, many nodes are in junior-default, but many nodes are in 
senior-default if and only if, loosely speaking, there are many more senior loans than junior ones. [Look at the dark 
region in the left-hand sides of Figs. [70) and 70).] In practice, we are most concerned with whether the number of 
bankruptcies of any type is large. Thus, our definition of the optimal seniority ratio in Sec. |II C] considered only the 
multiplex cascade region and not the fractions and 


Appendix E: Multilevel cascades under alternative settings 

In Appendix [C[ we claimed that the equality A max (J) = tr J, which is proved in the main text for our model, does 
not necessarily hold true in other variants of the multilevel cascade model. In what follows, we substantiate this claim 
by examining two variants of the 2-layer multiplex model. 


1. Exogenous senior thresholds 

In the model described in the main text, the threshold of senior-default, Rs, depends on the threshold of junior- 
default, Rj, and on the in-degree of the node on the junior layer, bj. In this subsection, we consider an alternative 
model of multilevel cascades in which the threshold of senior-default is an exogenous parameter. Let the (locally 
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fixed points (d> 4 ,, (b y J of the recursion equations 


(a) 


M 
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T3 



□ multiplex cascade: A max (J) > 1 

□ junior-only cascade: ^ jj > 1 

□ senior-only cascade: > 1 


0 2 4 6 8 



expected number of junior loans per bank = (/ j) 


expected number of junior loans per bank = (lj) 


numerical simulations 
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expected number of junior loans per bank = (lj) 



□ multiplex cascade: A max C7) > 1 

□ junior-only cascade: Jjj > 1 

□ senior-only cascade: > 1 


simulated (J) J 


(d) 
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□ multiplex cascade: A max (JJ) > 1 

□ junior-only cascade: Sfjj> 1 

□ senior-only cascade: Sss > 1 
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expected number of junior loans per bank = (lj) 


FIG. 7. (Color online) Here we compare the junior-only, senior-only and multiplex cascade conditions [Eqs. ([6j. 0, and ([sj, 
respectively] with the fixed points of the recursion equations 0 [in panels (a) and (b)] and with numerical simulations on 
random graphs with 10 4 nodes [in panels (c) and (d)[. The junior-default threshold of every node is Rj = 0.18 (as in Fig. 0. 
The axes are the edge densities (lj) and (Is) (i.e, the average numbers of junior and senior loans per bank). The solid blue, 
dashed orange, and dot-dashed green lines surround the multiplex, junior-only, and senior-only cascade regions [i.e., the sets 
of values of ((lj), (Is)) that satisfy the inequalities 0. 0, and ®, respectively]. In panels (a) and (b), the grayscale in the 
background shows the fixed points (j>io and (j^ of the recursion OT, begun from the initial condition (pg = (fig = 5 x 10 -4 . In 
panels (c) and (d), the grayscale in the background shows the fraction of nodes in junior- and senior-default, respectively, in 
numerical simulations with 10 4 nodes and 5 banks initially in senior-default, averaged over 75 simulations. In the upper-left 
part of all four panels, where most loans are senior, the fractions of nodes (4>ia,<t>oo) in junior- and senior-default at the end 
of the cascade are both « 1. By contrast, in the lower-right part, where most loans are junior, the fraction <p of nodes in 
junior-default is large [</>^, « 1; panels (a) and (c)[, whereas the fraction <j>%^ of nodes in senior-default is « 0.4 [see the light-gray 
region surrounded by the dashed orange line in panels (b) and (d)]. 
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treelike) graph hang down like a tree from a root node chosen uniformly at random. The recursion equation for the 
probability of default of a node t hops above the leaves of the tree is 


A“ — a( c 
h+i ~ 9 


>(#,#) = & + ( 1-&) £ 


J.out S',out 

Pi Pis 


E 


Is 


h+ls> 1 


rrij—O ms =0 


(El) 


for a £ {J, S'}, where p/) out and pf} out denote the out-degree distributions on the junior and senior layers, respectively. 
The response functions are now F a (l,m) — 1 mj+m s >R a (ij+i s ) f° r a ^ {^S}, where 0 < Rj < Rs < 1. 

The first-order cascade condition is thus given by A max ( l 7) > 1, where the entries of J are E[Zjli> ii .(; J+ ; g )] for 
(i,j) £ {J, S} 2 . Notice that the equality det J = 0 does not necessarily hold if Rj < Rs- It follows that A max (j) is 
not in general equal to trj 7 . 


2. Undirected edges 


Next, consider a multilevel cascade model in which edges are undirected. Again, let the graph hang down like a 
tree from a root node chosen uniformly at random. Consider a node located t hops above the leaves of the tree. 
The probability that this node is in a-default due to its children, conditioned on its parent in layer a £ {J, S} being 
solvent at layer a, is given by the recursion equations 


% + i =% + (!-%) E IjjPlY. E B l ^)B^ s (ef)Fj(l,m), 


lj-\-ls> 1 


h Is 

C E 

mj= 0 7715=0 


P?s l S 


flf +1 = 00 + (1 - Oq) E Phjfj E E B l ^)B l ^ s (ef)Fs(l,m), 

Zj+/S> 1 ' S' rn J=0 7715=0 


(E2a) 

(E2b) 


where one uses Eq. (E2aJ [respectively, Eq. (E2b)] if the node’s parent is a neighbor in layer J (respectively, layer S), 
and is the degree distribution on layer a. Because the graph is unweighted, and because this node lies at the end 
of an a-type edge chosen uniformly at random, its degree in layer a is the excess degree distribution, pf l a /{l a ). The 
response functions are 


Fj(l, fh) 
F s (l,fh) 


fl 

if > Rj 


= l 

(E3a) 

1 ° 

otherwise 

J 1 

if > Rj 

(E3b) 

|0 

otherwise 


Because the graphs are undirected, the term 
The corresponding Jacobian matrix is 


—bj in Fs in the original model [Eq. (2b)] is now —lj in Eq. (E3b). 


'v pi EsM l>i?j(Z J+Z5)] Po Mils-1)1 1 >Rj(lj+ls)\) 


where the expectations Ej and E 5 are over the joint probability distributions pf g and p J lj , respectively. It 
follows that det j" is generally nonzero, and therefore A max (jf) = tij does not necessarily hold. 


REFERENCES 


[1] P Gai and S Kapadia. Contagion in financial networks. Proceedings of the Royal Society A: Mathematical, Physical and 
Engineering Sciences, 466:2401-2423, June 2010. doi: 10.1098/rspa.2009.0410 















15 


[2] Erlend Nier, Jing Yang, Tanju Yorulmazer, and Amadeo Alentorn. Network models and financial stability. Journal of 
Economic Dynamics and Control, 31(6):2033-2060, 2007. URL: http://www.sciencedirect.com/science/article/pii/ 
SO165188907000097 doi:10.1016/j.jedc.2007.01.014 

[3] T. Kobayashi. A model of financial contagion with variable asset returns may be replaced with a simple threshold model 
of cascades. Economics Letters, 124:113-116, 2014. doi : 10.1016/j . econlet. 2014.05.003 

[4] P Gai, A Haldane, and S Kapadia. Complexity, concentration and contagion. Journal of Monetary Economics, 58:453-470, 

July 2011. doi:10.1016/j.jmoneco.2011.05.005 

[5] C. Upper. Simulation methods to assess the danger of contagion in interbank markets. Journal of Financial Stability, 
7:111 125, August 2011. doi: 10.1016/j . jfs.2010.12.001 

[6] R. M. May and N. Arinaminpathy. Systemic risk: the dynamics of model banking systems. Journal of the Royal Society 
Interface, 7:823-838, May 2010. doi: 10.1098/rsif . 2009.0359 

[7] T R Hurd and J.P. Gleeson. On Watts’ cascade model with random link weights. Journal of Complex Networks, l(l):25-43, 
May 2013. doi : 10.1093/comnet/cnt003 

[8] L. Bargigli, G. di Iasio, L. Infante, F. Lillo, and F. Pierobon. The multiplex structure of interbank networks. Quantitative 
Finance, 15(4):673-691, 2015. doi: 10.1080/14697688.2014.968356 

[9] M. Montagna and C. Kok. Multi-layered Interbank Model for Assessing Systemic Risk. Kiel Working Papers no. 1873, 
September 2013. 

[10] Matthew Elliott, Benjamin Golub, and Matthew O. Jackson. Financial networks and contagion. American Economic 
Review, 104(10):3115-53, 2014. URL: http://www. aeaweb. org/articles .php?doi=10.1257/aer. 104.10.3115 doi: 10. 
1257/aer.104.10.3115 

[11] Mikko Kivela, Alex Arenas, Marc Barthelemy, James P. Gleeson, Yamir Moreno, and Mason A. Porter. Multilayer networks. 
Journal of Complex Networks, 2(3):203-271, 2014. doi:10.1093/comnet/cnu016 

[12] S. Boccaletti, G. Bianconi, R. Criado, C.I. del Genio, J. Gomez-Gardenes, M. Romance, I. Sendina-Nadal, Z. Wang, 
and M. Zanin. The structure and dynamics of multilayer networks. Physics Reports, 544(1):1-122, 2014. URL: http: 
//www.sciencedirect.com/science/article/pii/S0370157314002105 doi:10.1016/j.physrep.2014.07.001 

[13] Sergey V Buldyrev, Roni Parshani, Gerald Paul, H Eugene Stanley, and Shlomo Havlin. Catastrophic cascade of failures 
in interdependent networks. Nature, 464(7291):1025-1028, April 2010. doi:10.1038/nature08932 

[14] Seung-Woo Son, Golnoosli Bizhani, Claire Christensen, Peter Grassberger, and Maya Paczuski. Percolation theory on 
interdependent networks based on epidemic spreading. EPL (Europhysics Letters), 97(1):16006, January 2012. doi: 
10.1209/0295-5075/97/16006 

[15] G J Baxter, S N Dorogovtsev, A V Goltsev, and J F F Mendes. Avalanche Collapse of Interdependent Networks. Physical 
Review Letters, 109(24):248701-5, December 2012. doi: 10.1103/PhysRevLett. 109.248701 

[16] Davide Cellai, Eduardo Lopez, Jie Zhou, James Gleeson, and Ginestra Bianconi. Percolation in multiplex networks 
with overlap. Physical Review E, 88:052811, Nov 2013. URL: http://link.aps.org/doi/10.1103/PhysRevE.88.052811 
doi:10.1103/PhysRevE.88.052811 

[17] Byungjoon Min, Su Do Yi, Kyu-Min Lee, and K.-I. Goh. Network robustness of multiplex networks with interlayer degree 
correlations. Physical Review E, 89:042811, Apr 2014. URL: http://link.aps.org/doi/10.1103/PhysRevE.89.042811 
doi : 10.1103/PhysRevE.89.042811 

[18] Charles D Brummitt, Kyu-Min Lee, and K-I Goh. Multiplexity-facilitated cascades in networks. Physical Review E, 
85:045102(R), April 2012. doi : 10.1103/PhysRevE. 85.045102 

[19] Osman Yagan and Virgil Gligor. Analysis of complex contagions in random multiplex networks. Physical Review E, 
86(3):036103, September 2012. doi: 10.1103/PhysRevE.86.036103 

[20] Kyu-Min Lee, Charles D Brummitt, and K-I Goh. Threshold cascades with response heterogeneity in multiplex networks. 
Physical Review E, 90(6), December 2014. doi: 10.1103/PhysRevE.90.062816 

[21] Duncan J. Watts. A simple model of global cascades on random networks. Proceedings of the National Academy of Sciences, 
99(9):5766-5771, 2002. doi: 10.1073/pnas .082090499 

[22] JP Gleeson and DJ Cahalane. Seed size strongly affects cascades on random networks. Physical Review E, 75(5):56103, 

2007. doi:10.1103/PhysRevE.75.056103 

[23] James P. Gleeson. Cascades on correlated and modular random networks. Physical Review E, 77(4):46117, 2008. doi: 
10.1103/PhysRevE.77.046117 

[24] Joshua Payne, Kameron Harris, and Peter Dodds. Exact solutions for social and biological contagion models on mixed 
directed and undirected, degree-correlated random networks. Physical Review E, 84(1), July 2011. doi : 10.1103/PhysRevE. 
84.016110 

[25] Sergey Melnik, Jonathan A Ward, James P Gleeson, and Mason A Porter. Multi-stage complex contagions. Chaos: An 
Interdisciplinary Journal of Nonlinear Science, 23(1):013124, 2013. doi: 10.1063/1.4790836 

[26] Michael Boss, Helmut Elsinger, Martin Summer, and Stefan Thurner. Network topology of the interbank market. Quan¬ 
titative Finance, 4(6):677-684, 2004. doi : 10.1080/14697680400020325 

[27] G. De Masi, G. Iori, and G. Caldarelli. Fitness model for the italian interbank money market. Phys. Rev. E, 74:066112, 
Dec 2006. URL: http://link.aps.org/doi/10.1103/PhysRevE.74.066112, doi : 10.1103/PhysRevE. 74.066112 

[28] Rama Cont, Amal Moussa, and Edson B Santos. Network structure and systemic risk in banking systems. In Jean-Pierre 
Fouque and Joseph A. Langsam, editors, Handbook on Systemic Risk. Cambridge University Press, New York, 2013. 

[29] L C G Rogers and LAM Veraart. Failure and Rescue in an Interbank Network. Management Science, 59(4):882-898, 
April 2013. doi : 10.1287/mnsc . 1120.1569 

[30] Helmut Elsinger. Financial Networks, Cross Holdings, and Limited Liability . Oesterreichische Nationalbank Working 




16 


Paper 156, May 2009. 

[31] C Gourieroux, J C Heam, and A Monfort. Liquidation equilibrium with seniority and hidden CDO. Journal of Banking 
and Finance, 37(12):5261-5274, December 2013. doi: 10.1016/j . jbankfin.2013.04.016 

[32] L. Eisenberg and T. Noe. Systemic risk in financial systems. Management Science, 47(2):236-249, 2001. doi: 10.1287/ 
mnsc.47.2.236.9835 

[33] D Centola, V Egufluz, and M W Macy. Cascade dynamics of complex propagation. Physica A: Statistical Mechanics and 
its Applications, 374(l):449-456, January 2007. doi:doi: 10.1016/j .physa. 2006.06.018 

[34] Joshua Payne, Peter Dodds, and Margaret Eppstein. Information cascades on degree-correlated random networks. Physical 
Review E, 80(2), August 2009. doi: 10.1103/PhysRevE.80.026125 

[35] R R Liu, W X Wang, Y C Lai, and B H Wang. Cascading dynamics on random networks: Crossover in phase transition. 
Physical Review E, (85):026110, 2012. doi: 10.1103/PhysRevE.85.026110 

[36] Mark Granovetter. Threshold models of collective behavior. American Journal of Sociology, 83(6): 1420-1443, 1978. URL: 
http://www.jstor.org/stable/2778111 

[37] G. Bianconi. Statistical Mechanics of Multiplex networks: Entropy and Overlap. Physical Review E, 87:062806, 2013. 
doi:10.1103/PhysRevE.87.062806 

[38] Data on leverage ratios of banks was accessed at http://www.bankregdata.com/allHMmet.asp?met=LEV on January 22, 
2015. The source of the data is the Federal Financial Institutions Examination Council (FFIEC), available at https: 
//cdr.ffiec.gov/public/ 

[39] K-I Goh, B Kahng, and D Kim. Universal Behavior of Load Distribution in Scale-Free Networks. Physical Review Letters, 
87(27):278701, December 2001. doi: 10.1103/PhysRevLett.87.278701 

[40] M Catanzaro and R Pastor-Satorras. Analytic solution of a static scale-free network model. The European Physical Journal 
B, 44(2):241-248, April 2005. doi:10.1140/epjb/e2005-00120-9 

[41] The original static model introduces correlations in the degrees of neighbors j40]. To avoid such correlations, we constrain 
degrees to be less than y/N (where N is the number of banks) [JS], and we use the configuration model with in- and 
out-degrees drawn from the analytical solution of the degree distribution in the static model 1401 Eq. (35)]. 

[42] Nicholas Beale, David G. Rand, Heather Battey, Karen Croxson, Robert M. May, and Martin A. Nowak. Individual versus 
systemic risk and the regulator’s dilemma. Proceedings of the National Academy of Sciences, 108(31): 12647-12652, 2011. 
doi:10.1073/pnas.1105882108 

[43] Charles D Brummitt, Rajiv Sethi, and Duncan J Watts. Inside Money, Procyclical Leverage, and Banking Catastrophes. 
PloS one, 9(8):el04219, August 2014. doi: 10.1371/journal .pone. 0104219 

[44] Fabio Caccioli, J. Doyne Farmer, Nick Foti, and Daniel Rockmore. Overlapping portfolios, contagion, and financial stability. 
Journal of Economic Dynamics and Control, 51(0):50-63, 2015. URL: http://www.sciencedirect.com/science/article/ 
pii/SO165188914002632, doi:10.1016/j.jedc.2014.09.041 

[45] Michele Catanzaro, Marian Boguna, and Romualdo Pastor-Satorras. Generation of uncorrelated random scale-free net¬ 
works. Physical Review E, 71(2):027103, February 2005. doi: 10.1103/PhysRevE.71.027103 




